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Zhi-zhong Xing %° Zhen-hua Zhao ° * 


“Institute of High Energy Physics & School of Physical Sciences, University 
of Chinese Academy of Sciences, Beijing 100049, China 
’Department of Physics, Liaoning Normal University, Dalian 116029, China 
“Center for High Energy Physics, Peking University, Beijing 100080, China 


Abstract 


The neutrinoless double-beta (0v2) decay is currently the only feasible process in particle and 
nuclear physics to probe whether massive neutrinos are the Majorana fermions. If they are of the 
Majorana nature and have a normal mass ordering, the effective neutrino mass term (mM) se of a 0v28 
decay may suffer significant cancellations among its three components and thus sink into a decline, 
resulting in a “well” in the three-dimensional graph of |(m),,| against the smallest neutrino mass 
m, and the relevant Majorana phase p. We present a new and complete analytical understanding 
of the fine issues inside such a well, and identify a novel threshold of |(m),,| in terms of the neutrino 
masses and flavor mixing angles: |(m),.|, = Mmg sin? 0,3 in connection with tan 043 = \/m,/m, and 
p=. This threshold point, which links the local minimum and maximum of |(m),.|, can be used 
to signify observability or sensitivity of the future 01v26-decay experiments. Given current neutrino 


oscillation data, the possibility of |(m)..| < (m) is found to be very small. 
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Since Ettore Majorana first formulated a fermionic particle that should be its own antiparticle in 
1937 [1], a huge amount of attention has been paid to the Majorana fermions in particle and nuclear 
physics and the Majorana zero modes in solid-state physics [2]. In particular after the experimental 
discoveries of solar, atmospheric, reactor and accelerator neutrino oscillations [8], whether massive 
neutrinos are the Majorana fermions becomes an especially burning question among a number of 
fundamentally important questions in neutrino physics and cosmology. If this is the case, then the 
neutrinoless double-beta (072) decays of some even-even nuclei are expected to take place [4]. Namely, 
N(A, Z) > N(A, Z +2) + 2e7, where the lepton number is violated by two units. Given the fact that 
the neutrino masses are so small that all the lepton-number-violating processes must be desperately 
suppressed, currently the unique and only feasible way to demonstrate the Majorana nature of massive 
neutrinos is to observe the 0v26 decays. In this respect a number of ambitious experiments are either 
underway or in preparation [5]. 

In the standard scheme of three neutrino flavors the rate of a 0v2( decay is proportional to the 
squared modulus of the effective Majorana neutrino mass term [6] 


(M) ce = Mm [Ua e? + mUe]? + m3|Ueg|"e'7 ’ (1) 


where m; denotes the i-th neutrino mass (for i = 1,2,3), U,; is the corresponding element of the 
3 x 3 neutrino mixing matrix U [8], and p and ø stand for the Majorana phases. One often chooses 
to parametrize |U,,| as follows [3]: |U.,| = cos 015 cos 6,3, |U.9| = sin 6,9 cos 6,3, and |U,3| = sin 6,5. 
The three mixing angles 6,5, 0,3 and əz have been determined to a good degree of accuracy from 
current neutrino oscillation data, so have been the value of Am3, = m3 — m? and the modulus of 
Am3, = m} — mî [B]. But the sign of Am3, and the two phase parameters in Eq. (1) remain unknown, 
nor does the absolute neutrino mass scale. That is why |(m),,| is usually plotted as a function of m, 
in the normal mass ordering (NMO) case (Am3, > 0) or ms in the inverted mass ordering (IMO) case 
(Am3, < 0) by allowing p and o to vary from 0 to 27 [9]. In such a so-called Vissani graph, a two- 
dimensional “well” can appear in the NMO situation due to a significant cancellation among the three 
components of (m)... The bottom of the well signifies the case of |(m),.| — 0 [10], a disappointing 
possibility which is definitely consistent with the present experimental data. 

Two immediate questions are in order: (1) how possible for the three neutrinos to have a NMO; (2) 
how possible for the actual value of |(m)_..| to fall into the well and become unobservable in any realistic 
0v28 experiments. A combination of current atmospheric (Super-Kamiokande [11]) and accelerator- 
based (T2K and NOvA [13]) neutrino oscillation data preliminarily favors the NMO at the 2c 
level. If this turns out to be the case, an answer to the second question will be highly desirable because 
it can help interpret the discovery or null result of a 0v26 experiment in the standard three-flavor 
scheme, although some kind of hypothetical (ad hoc) new physics may also contribute to |(m),.|. 

The present work aims to answer the second question by giving a new and complete analytical 
understanding of the fine structure of the three-dimensional well of |(m),,| against m, and p, as 
illustrated in Fig. 1, where the best-fit values Am3?, = 7.54 x 10~° eV’, Am, = 2.47 x 1073 eV?, 
sin? 0,5 = 0.308 and sin? 013 = 0.0234 have been taken as the typical inputs. We identify a novel 
threshold of |(m),.| which is located at the center of the well: |(m),.|, = Mmg sin? 0,3 in connection with 
tanĝis = \/m,/m, and p = r. This threshold point links the local minimum and maximum of | (m) 


and it can be used to signify the observability or sensitivity of the future 0v2G-decay experiments. 


eal 


Given current neutrino oscillation data, the possibility of |(m)..| < |(m)..|, is found to be very small. 


bale 


'The phase convention taken here is highly advantageous when considering the interesting and experimentally-allowed 
neutrino mass limit mı > 0 (or m3 —> 0), in which p (or o) automatically disappears [7]. 
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as functions of m, and p in the NMO case, where the best-fit values Ams, = 7.54 x 1075 eV?, 


Am, = 2.47 x 107? eV’, sin? 6,5 = 0.308 and sin? 6,, = 0.0234 have typically been input. 


Figure 1: Three-dimensional illustration of the upper (orange) and lower (blue) bounds of |(m) 


Fig. 1 shows that the depth of the well of |(m),..| is mainly sensitive to a narrow parameter space 
of m, and p, while the other Majorana phase o plays an important role in shaping the bottom of the 


o 
well [I5]. The latter point can be seen in an analytical way as follows. Taking Fel = 0, we obtain 
loj 


mı sin p 


(2) 


tano = 


m, cos p + Mo tan? bjs ’ 
so as to maximize or minimize |(m}ee| for the given values of m, and p. Substituting Eq. (2) into the 
expression of |(m)}ee| in Eq. (1), one arrives at the following upper (“U”) and lower (“L”) bounds: 


Km) eelu = [m2 cos? 013 + mg sin? O13] (3) 


where the sign “+” (or “—”) corresponds to “U” (or “L”), and 


1 
My = yr cos* 6,5 + zma sin? 2043 cos p + MŽ sinf 6,5 . (4) 


It is easy to understand this result in an intuitive way: for any given values of m, and p, the maximum 
of |(M)ee| comes out when the sum of the first two components of (m),, has the same phase as the 
third one (i.e., o); and the minimum of |(m),.| arises when the difference between these two phases is 


equal to +r. The bottom of the well shown in Fig. 1 corresponds to |(m),.|;, = 0, or equivalently 


Given the expressions my = 4/ M? + Am3, and mz = \/m? + Am, in the NMO case, Eq. (5) allows 
2 1 21 3 1 31 


us to fix how the two free parameters m, and p are correlated with each other. Using the same 
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Figure 2: The numerical correlation between m, and p in three typical cases: (a) |(Mm)eelp = 0 


(the red curve); (b) Km)eelg = (Meel = Mz sin? 6,3 (the black dot and the blue curve); and (c) 
Kmecl = n|(Meelą with n > 2 (the black curves). Here the best-fit values of Am3,, Am2}, sin? 643 


and sin? 6,3 used in plotting Fig. 1 have been input. 


best-fit inputs of Am3,, Am3,, sin? 6,5 and sin? 6,3 as those used in plotting Fig. 1, we illustrate 
the numerical correlation between m, and p dictated by Eq. (5) in Fig. 2 — the red curve. Such a 
correlation curve roughly looks like an ellipse, but a careful analytical check shows that it does not 
really obey the standard equation of an ellipse. Fig. 2 tells us that touching the bottom of the well 
(i.e., |(m),..| + 0) is not a highly probable event at all, because it requires m, and p to lie in the 
narrow regions 2 meV £ m; <7 meV and 0.86 < p/m < 1.14, respectively [16]. 

Another salient feature of the well is the “bullet”-like structure of |(m),.|;, as shown in Fig. 1, 
corresponding to the parameter space of Mis < Ma tan? 6,3. In other words, the surface of this bullet 
is described by 


(Mm) eel, = ™3 sin? 013 — Mi cos? CER (6) 


The extremum of |(m),.|;, in this inner region of the well is supposed to be located at a point fixed 


by the following two conditions: 


o |M) eel, Mmo sin? 2943 cos? 3. 
= si 


np=0 
Op AM 5 p : 
2 2 o2 D 2 
O\(M) eel, Mi . 2 MoM? — AM3; sin? 642 (mg sin? 04a + M4 cos? Ojo cosp) v E 
= sin“ 013 — cos“ 643 = 0. (7) 
dm, m3 Meia Meiz 


The first condition definitely leads us to p = 0 or m. But Fig. 2 clearly shows that p should only take 
a value around ~r inside the well, and thus it is appropriate to take p = ~m instead of p = 0. In this 
case Mijo = |m] cos? 045 — Mo sin? 043| holds, and the second condition in Eq. (7) is simplified to 


O (Meelu _ Mı -2 2 as 
= sin? 0,3, + | cos? 04o — — sin? 6,4 | cos? 6,, = 0 8 
dm, Mg 13 2 m 12 13 ’ (8) 
where “+” correspond to the prerequisites m4 < my tan? 64a and m, > mz tan? 0,5, respectively. But 


in reality Eq. (8) can never be fulfilled since its second term is much larger than its first term as a 
result of (a) 2.50 x 1071 < sin? 04 < 3.54 x 107! and 1.85 x 107? < sin? 043 < 2.46 x 107? at the 30 


level and (b) m,/m3 < m,/mp in the NMO case. Nevertheless, Eq. (8) can at least allow us to 
draw a conclusion that is absolutely consistent with current experimental data: 


ð 
Ə |m elr, 2 0 for mM < Mo tan? O19 5 
my 
ð | (M) eel 
Z a <0 for m > mtan? bz. (9) 


This observation means that |(m),,|;, increases when m} < my tan? 6; holds, and it decreases when 


mı > mM, tan? 6,5 holds. Hence there must be a local maximum for |(m),.|;,, denoted as 


(7) cel = mg sin? 613 = Vy mi T Am; sin? 013 (10) 


at the position fixed by p = 7 and 


23 
2 l2 2 2 [n2 S O12 
Mı = M tan" O49 = 4/ MÍ + Ams, tanĝi => mı = 4/AmMm5 Joos Oy : (11) 
cos 2645 


In Fig. 1 this point is exactly the tip of the bullet inside the well! In other words, the local maximum 
of (Meel, arises from Eq. (6) at Mio = 0. Given the best-fit values of Am3,, Am3,, sin? @,. and 
sin? 0,; that have been used in plotting Fig. 1, the numerical location of the tip of the bullet turns 
out to be (m4, p, |(M})eelx) =~ (4 meV, 180°, 1 meV). 

The above analysis explains why the bottom of the well does not converge to a single point and why 
it is not flat either. In a similar way one can understand why there is a local minimum for |(m),¢|y, 


as shown in Fig. 1. The extremum of |(m),,|y is expected to be located at a position determined by 


I|(M)cely _ Mym sin? 20,5 cos? 643 


sinp=0, 
Op Amy 5 P 


IA\(M) eely — M sin? 013+ = 
Om, Ms MMM 


MMi — Am’, sin? 843 (M sin?” 642 + M4 cos” 6,5 cos p) cos? 8,4 = 0. (12) 
ig ss 


Of course, only p = 7 is allowed with respect to the first condition in Eq. (12). The second condition 
in Eq. (12) can never be satisfied for the same realistic reasons given below Eq. (8). An analogous and 
straightforward analysis tells us that the local minimum of |(m),,|y exactly coincides with the local 
maximum of |(m),.|;,, and thus both of them are described by Eqs. (10) and (11). This interesting 
result explains why the upper (in orange) and lower (in blue) bounds of |(m),,| connect with each 
other in Fig. 1 when m, = ym? + Am3, tan? 0j and p = m hold. Note that the overlap of the local 
maximum of |(m),.|;, and the local minimum of |(m),.|y can also be understood from Eq. (3) itself. 
At my = ym? + Am, tan? 6, and p = 7, one simply has |(m),.|;, = |(m)eely = mM3sin? 0; as a 
consequence of Mja = 0. So |(m)..|, = m3 sin? 043 ~ 1 meV stands for a threshold of |(m),,| in the 
NMO case. 

To visualize the steepness of the slope of |(m),,.|;, around the well in Fig. 1, let us project its 
contour onto the m,-p plane by taking |(m),.|;, = n|(m),-|, (for n = 0,1,2,---) in Fig. 2. It is 
especially interesting to compare between the contours of the well at its bottom with |(m),.|;, = 0 
(the red curve) and at its threshold height with |(m)..|;, = |(m) the blue curve and the black 
point). They clearly show how the well becomes narrower when the value of |(m),.|;, goes down. The 


ls ( 


profile of |(m),.|;, will be partially open and thus lose its “well” feature as | (meel > 2|(m) eel, is taken 


zely 
into account. Now that |(m)eelg > |(™m)¢el, always holds outside the blue curve in Fig. 2, we argue 
that the parameter space of |(mMm)eelg < |(M)ce|, (i-e., 0.4 meV Sm, < 10 meV and 0.66 S p/m < 1.34) 


is a simple measure of the chance for |(m),.| to fall into the well and become completely unobservable. 
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Figure 3: The parameter space of m}, p and o allowed for |(m),..| < |(M)eelą to hold, where the best-fit 
values of Am3,, Am?3,, sin? 6, and sin? 0,3 used in plotting Fig. 1 have been input. The intersecting 
surfaces for m, = 1,2,4 and 6 meV on the p-o plane are explicitly shown in the figure. 


In general, |(m),.| depends on all the three unknown parameters m}, p and ø. To illustrate how 
probable or improbable for |(m),..| to have a value smaller than |(m),,|, in a more explicit way, we plot 
the three-dimensional parameter space of m,, p and ø in Fig. 3, where the best-fit values of Am3,, 
Am2, sin? 0, and sin? 0}; used in plotting Figs. 1 and 2 have been input. For clarity, the intersecting 
surfaces on the p-o plane corresponding to m, = 1, 2,4 and 6 meV are specified in the figure. One can 
see that this parameter space is very small as compared with the whole cubic space (i.e., the whole 
regions of m,, p and o allowed by current experimental constraints). In comparison with m, and p, 
the phase o is only weakly constrained in Fig. 3. When the first two components of (m)ee in Eq. (1) 
essentially cancel each other out (i.e., 2 meV Sm, <7 meV and 0.86 < p/m < 1.14), a large part of 
the range of o is allowed (e.g., the black intersecting surface corresponding to m, = 4 meV in Fig. 3). 
But when the value of m, decreases, the value of ø should approach 7, such as the green intersecting 
surface corresponding to mı = 1 meV in Fig. 3. In this case the second component of (m),,. in 
Eq. (1) can be cancelled by the other two components to a maximal level. For a similar reason, 
the value of ø should approach 0 or 27 when the value of m} increases (e.g., the blue intersecting 
surface corresponding to m, = 6 meV in Fig. 3). In any case we conclude that the possibility of 
|(m)ee| < |(m),,-|, involves significant cancellations among its three components and is really small. 

From an experimental point of view, the threshold |(m),.|, should signify an ultimate limit of 
the reachable sensitivity to |(m)..| in the foreseeable future. At present the most sensitive 0v2(- 
decay experiments can only set an upper limit of |(m),,.| around 165 meV [17], which depends on 
some theoretical uncertainties in calculating the relevant nuclear matrix elements [18]. The most 


ambitious next-generation high-sensitivity 0v23-decay experiments (e.g., nEXO [19]) are likely to 


probe |(m}ee| at the level of a few tens of meV P| [5], a sensitivity still much larger than the threshold 
value |(m)..|, ~ 1 meV ]?| In this sense there would be no hope to observe any 0v2(6-decay signal if 


in the standard three-flavor scheme. 


|(m).-| were unfortunately around or below the value of |(m)..|, 


Before ending our discussions about (m),. and its possible parameter space in the NMO case, let 
us briefly comment on the relationship tan 6,5 = y m/m from a model-building point of view. This 
condition, together with p = 7, allows for |(m),.| = KM) eels = Mg sin? 6,3 as a remarkable threshold. 
It is well known that the Cabibbo angle @, of quark flavor mixing can be related to the ratio of 
quark masses mg and m, in a class of models [22]: tan 6c ~ VMg Ms , which is consistent with the 
experimental data to a good degree of accuracy. In comparison, the possibility of tan 6,5 ~ ym /m 
is also interesting, in particular when the NMO is true for the three mass eigenstates of vs, Vy and v 
neutrinos. For example, we find that an effective Majorana neutrino mass matrix of the form 


0 AA a V2 sinOj, +i —i 
M,=|A BC -ma +i 0 ol, (13) 
ACB =i 0 0 


where A, B and C are all real, can essentially predict |(m),.| = Mg sin? 6,3 and tan@,. = \/m,/m, 
together with 6,3 = 7/4, 6 = —7/2, p = 7 and o = 0 in the standard parametrization of U. Because 
M, possesses the exact u-T reflection symmetry, which can easily be simplified to the u-r permutation 
symmetry in the 6,3; — 0 limit, one may take it as a starting point to build a phenomenological 
neutrino mass model in this connection 

In summary, we have achieved some new and important insights into the effective neutrino mass 
(M)ee of the 0v28 decays in the NMO case — a case which seems to be more likely than the IMO 
case according to today’s preliminary experimental data. Because |(m),.,| depends not only on the 
unknown neutrino mass m but also on the free Majorana phases p and øg, a novel three-dimensional 
presentation of |(m),.| against m, and p reveals an intriguing “well” structure in the NMO case. 
The present work provides a new and complete analytical understanding of the fine issues inside such 
a well. We find a particularly interesting threshold of |(m)..| in terms of the neutrino masses and 
|e = M3 sin? #,3 in connection with tan @,. = Jfm,/m and p = 7. We 
suggest that this threshold point, which links the local minimum and maximum of |(m)..|, be used to 


flavor mixing angles: |(m),. 
signify observability or sensitivity of the future 0v26-decay experiments. In view of current neutrino 


must be very small. In other 


oscillation data, we conclude that the possibility of |(m)..| < KM)eelx 


words, it should be very promising to detect a signal of the 0726 decays and verify the Majorana 
nature of massive neutrinos in a foreseeable future, even if they have a normal mass spectrum. 
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2Note that the accuracy of a prediction for the experimental sensitivity crucially depends on our knowledge of the 
relevant nuclear physics. In the worst possible scenario, uncertainties from nuclear physics might even weaken the 
expected experimental sensitivities by a factor as large as 5 [5]. 

3In Ref. [20] a purely statistical analysis of the possibility of |(m)..| <1 meV has been done to see to what extent the 
Majorana phases p and o can be constrained for a given value of mı. While in Ref. the conditions for |(mMm)ee| > 1 
meV are analyzed in the special case of mı — 0 or 043 > 0. 
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